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Abstract 

We present results on the production of high transverse momentum charm 
mesons in collisions of 515 GeV/c negative pions with beryllium and copper 
targets. The experiment recorded a large sample of events containing high 
transverse momentum (pr) showers detected in an electromagnetic calorime- 
ter. From these data, a sample of charm mesons has been reconstructed via 
their decay into the fully charged Kirn mode. A measurement of the single 
inclusive transverse momentum distribution of charged D mesons from 1 to 
8 GeV/c is presented. An extrapolation of the measured differential cross 
section yields an integrated D^ cross section of 11.4±2.7(sioi) ±3.3(syst) /ib 
per nucleon for D mesons with xf > 0. The data are compared with expec- 
tations based upon next-to-leading order perturbative QCD, as well as with 
results from pythia. We also compare our integrated D cross section with 
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measurements from other experiments. 

PACS number(s): 13.85.Ni, 14.40.Lb, 13.60.Hb, 13.85.-t 
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I. INTRODUCTION 

Over the past decade, measurements of charm production in hadronic interactions have 
provided valuable tests of the applicability of perturbative QCD (pQCD) to the production 
of heavy quarks [0-§]. Quark- ant iquark annihilation and gluon fusion are the leading order 
(LO) contributors to charm quark hadroproduction. The next-to-leading order (NLO) con- 
tributions to the cross sections have been evaluated |)|-|| , and are comparable to the leading- 
order contributions. While the results of the NLO calculations can accommodate the charm 
cross sections observed in 7r~-nucleon interactions p|-[T5|, the large size of the NLO contribu- 



tions is an indication that still higher order contributions may be significant. Furthermore, 
NLO calculations of the total charm quark cross section exhibit significant sensitivity to the 
choice of input parameters, including the charm quark mass (m c ), the renormalization and 
factorization scales, as well as the parton distribution functions (PDF's). For example, the 
calculated charm cross section changes by a factor of ~3 when the renormalization scale is 
varied from 2m c to \m c . Varying the charm quark mass from 1.2 GeV/c 2 to 1.8 GeV/c 2 
changes the calculated charm cross section by as much as an order of magnitude fl6| , |l7l . 
While there is significant theoretical uncertainty in the total charm quark production cross 
section, there is less uncertainty in the shapes of the differential distributions. The shapes 
of the LO and NLO single inclusive charm quark distributions versus px are rather similar, 
and their shapes exhibit smaller sensitivity to variations in m c or the renormalization scale 

IB 



To compare measurements of charm hadron production to the results of pQCD calcu- 
lations, the consequences of hadronization of the produced charm quarks must be taken 
into account. The fragmentation of charm quarks into charm hadrons is inherently a low 
momentum transfer process, and is therefore currently beyond the domain of pQCD. Nev- 
ertheless, the effects of fragmentation may be described phenomenologically by convoluting 
the partonic cross sections with a suitable fragmentation function. One hopes to describe the 
hadronization of charm quarks via a universal, pro cess- independent fragmentation function, 



such as the Peterson et al. |19| form, as measured in e + e~ collisions. Convoluting the NLO 
prediction for charm quark production with a fragmentation function results in a softening of 
the predicted p? spectrum of the charm hadrons relative to the charm quarks. The Peterson 
et al. fragmented NLO result is softer than the measured px spectrum for charm hadrons, 
and it has been observed that the unfragmented NLO result for charm quarks reproduces the 
shape of observed charm hadron pt spectra reasonably well in the kinematic range xp > 
and pt < 4 GeV/c [0J. One might expect additional nonperturbative effects, such as the 
intrinsic transverse momentum of the incoming partons, to have an impact on the measured 
differential distributions. Frixione et al. |TJ noted that the Peterson et al. fragmented NLO 
calculations of the charm pt spectra can be brought into agreement with data from exper- 
iments E769 and WA82, provided the partons (in each hadron) are supplemented with an 
average squared intrinsic transverse momentum, (k 2 r )=2.0 (GeV/c) 2 . 

It is also of interest to compare the data with a Monte Carlo (MC) simulation that 
incorporates a model of the fragmentation process (such as the LUND string model |2(J 
as implemented in pythia). Again, nonperturbative effects may need to be taken into ac- 
count to match the experimental results, pythia simulates the so-called "leading particle 
effect" , which results in an enhanced forward production for D mesons whose light quark is 



a spectator valence quark from the incoming beam particle. In addition to initial state radi- 
ation effects, pythia also includes additional intrinsic transverse momentum characterized 



by the parameter kj_ [2T|] . Reasonable agreement has been achieved between a pythia MC 
simulation and the E769 data using a kj_ value of ~1 GeV/c ||22|| . 

Fermilab experiment E706 |23j was designed to study high pt phenomena, principally 



associated with direct photons |24],[^J and high pt jets ||26|| . Jets arising from large mo- 
mentum transfer collisions are expected to be rich in high pt charm particles. The high pt 
requirement of the E706 trigger enhances the fraction of selected events containing charm by 
nearly an order of magnitude compared to a minimum bias trigger. Since the selected events 
result from high transverse momentum interactions, the data constitute a unique sample in 
which to study charm particles (other recent fixed target charm experiments frequently em- 
ployed lower threshold or minimum bias triggers, and yielded a rich sample of mostly lower 
Pt events). From a theoretical standpoint, one might expect the pQCD calculations to be- 
come more reliable in the kinematic range accessible to E706. This paper presents results 
of a study of high p? charm particles produced in 515 GeV/c 7r~-nucleon collisions. The 
measured differential cross sections are compared to NLO pQCD calculations and results 
from a pythia MC simulation. We also compare the integrated D ± cross section with NLO 
pQCD calculations and other recent measurements. 

II. APPARATUS 

Experiment E706 was performed in the Fermilab Meson West beam line. The unsepa- 
rated negative secondary 515 GeV/c beam was primarily composed of pions with a small 
admixture of kaons (<5%). Figure [l] displays a diagram of the key elements of the Meson 
West spectrometer for this measurement |[27|| . The detector included a precision charged par- 
ticle tracking system and a large acceptance liquid argon calorimeter (LAC). The charged 
tracking system employed silicon microstrip detectors (SSD), a large aperture dipole analysis 
magnet, proportional wire chambers (PWC), and straw tube drift chambers (STDC). The 
LAC contained a finely segmented electromagnetic section (EMLAC) as well as a hadronic 
section. The high p? trigger was based upon signals originating from showers detected in 
the EMLAC. Only data from the charged particle tracking system and the electromagnetic 
section of the calorimeter contributed directly to this analysis. 

A. Charged particle tracking system 

The target region of the Meson West spectrometer (shown in Fig. |2|) consisted of nuclear 
targets and a SSD system p8| . The targets included two 780 /im thick copper pieces followed 
by two beryllium cylinders of length 3.71 cm and 1.12 cm, respectively. A 1 cm air gap 
between the two beryllium cylinders was designed for use in searches for heavy quark decays 
|2l| . The targets were supported in a Rohacell stand which had a cylindrical hole bored 
along the beam axis where the targets were positioned. Three beam SSD modules were 
located upstream of the target and 5 vertex SSD modules were located downstream of the 
target. Each module was composed of a pair of single-sided SSD planes with strips aligned 
vertically and horizontally along the X and Y axes, respectively. All of the SSD planes 



were 50 //m pitch detectors, with the exception of the pair of vertex SSD's nearest the 
target, which consisted of a high resolution 25 /im pitch central region and 50 /zm pitch 
outer regions. A total of 8912 strips were instrumented, providing an angular acceptance of 
m ±125 mrad in each view. 

A large dipole analysis magnet was located downstream of the SSD system. Charged 
particles passing through the magnet received a pt impulse of ^450 MeV/c in the horizontal 
plane. Four PWC stations [§3'0] were located downstream of this dipole analysis magnet, 



each separated by ~1 m. Each station consisted of four proportional wire sense planes with 
wires oriented at the angles —90° (X view), 0° (Y view), 37° (U view), and —53° (V view). 
Thus, the X and Y views were orthogonal to one another, as were the U and V views. 
Within each sense wire plane, the spacing between the sense wires was 0.254 cm. A total 



of 13,440 sense wires were instrumented. The STDC system [32] consisted of two stations, 
with each station consisting of 4 X view planes followed by 4 Y view planes. The X and Y 
view planes in the upstream station consisted of 160 and 128 tubes, respectively. Each tube 
had a diameter of 1.03 cm. Each plane of the downstream STDC consisted of 168 tubes, and 
each tube had a diameter of 1.59 cm. The hit resolution of individual tubes in the STDC's 
was typically in the range of 200 to 300 /im. 

B. Electromagnetic calorimeter 

The EMLAC [p3lp4| , which was located 9 m downstream of the target, was assembled 
from four independent quadrants, each of which instrumented j of the azimuthal acceptance 
of the detector. The inner and outer radii of the EMLAC were 20 cm and 160 cm respectively, 
corresponding to a center of mass rapidity coverage of —1 to 1 (for the incident 515 GeV/c 
beam). Each quadrant was composed of 66 layers. Each layer consisted of an absorber 
sheet (which was 2 mm of lead in all but the first layer), a 2.5 mm liquid argon gap, a 
pair of octant-size copper-clad G-10 anode boards, and another 2.5 mm argon gap. In 
alternating layers, the copper cladding of the G-10 anode boards was cut to form either 
concentric (R) strips or azimuthal (0) strips. The locations and widths of the R strips on 
the radial anode boards were such that the R strips were focussed on the target region so that 
neutral particles produced in the target passed through the same R strip in each successive 
radial anode board. The width of the R strips on the first anode board was 5.5 mm. The 
interleaved azimuthal anode boards were subdivided at a radius of 40 cm into inner and 
outer regions. Each of the inner strips subtended an angle of ^ radians in azimuth 
while the outer strips subtended an angle of ^ radians. Longitudinally, the EMLAC was 
read out in two sections. The front section consisted of the first 22 layers (8.5 radiation 
lengths) while the back consisted of the remaining 44 layers (18 radiation lengths). For each 
section in each octant, signals from corresponding R (0) strips were ganged together and 
read out independently. 

C. Trigger 

The trigger selected events producing high transverse momentum showers in the EM- 
LAC. This event selection process involved four stages; beam and interaction definitions, a 



pretrigger requirement, and the final trigger requirements | 35f . The beam definition required 



that a single beam particle was detected in the beam hodoscope located 2 m upstream of 
the target. A scintillation counter with a | inch diameter hole was located downstream of 
the beam hodoscope and was used to reject interactions initiated by particles in the beam 
halo. Pairs of interaction counters were mounted near the upstream and downstream mirror 
plates of the dipole analysis magnet. An interaction was defined as a coincidence between 
signals from at least two of these four interaction counters. To minimize potential confu- 
sion due to out-of-time interactions, a cleaning filter rejected any interactions that occurred 
within ±60 ns of each other. For those interactions that satisfied the beam and interaction 
definition, the px deposited in various regions of the EMLAC was evaluated by weighting 
the energy signals from the EMLAC R channel amplifier fast outputs by ~ sin#j, where 9i 
is the polar angle that the i th strip subtends with respect to the nominal beam axis. The 
pretrigger p? requirement was satisfied if the px detected in the inner 128 R channels or 
the outer 128 R channels of at least one octant was greater than a threshold of ~1.7 GeV/c 
(for the PRETRIGGER HI). A pretrigger signal was issued only if the signals from a given 
octant satisfied that pretrigger pt requirement and there was no evidence in that octant of 
substantial noise or significant pt attributable to an earlier interaction and there was no 
incident beam halo muon detected [KG]. The pretrigger signal latched the data from the 
various subsystems while the final trigger decision was being evaluated. 

The experiment employed several different high p? trigger definitions that were based 
upon the LOCAL and GLOBAL signals from octants that satisfied the pretrigger. Local 
trigger groups were formed by clustering the 256 R channels in each octant into 32 groups 
of 8 channels. Each of the adjacent pairs of groups of 8 channels (groups 1&2, 2&3, . . . , 
31&32) formed a local group of 16 strips. If the pt detected in any of these groups of 16 was 
above a specified high (or low) threshold, then a LOCAL HI (or LO) signal was generated for 
that octant. A GLOBAL HI (or LO) signal for a given octant was generated if the summed 
Pt from the groups of 8 in that octant was above a specified high (or low) threshold. In 
order to suppress coherent noise effects, only groups of 8 registering at least ^250 MeV/c 
contributed to this global px sum. (This cutoff was applied independently to signals from 
groups of 8 from the front and back sections of the EMLAC.) Three of the trigger types, 
which accounted for ^^80% of the E706 data, were used in this charm analysis; they were 
the SINGLE LOCAL HI, the LOCAL®GLOBAL HI, and the TWO-GAMMA triggers. The 
SINGLE LOCAL HI trigger required a LOCAL HI signal from an octant that satisfied the 
PRETRIGGER HI. The LOCAL HI threshold was ^3 GeV/c. The LOCAL®GLOBAL HI 
trigger required the coincidence of a LOCAL LO signal (threshold «1.7 GeV/c), a GLOBAL 
HI signal (threshold ^3 GeV/c), and a PRETRIGGER HI all from the same octant fJ7|. 
The TWO-GAMMA trigger required LOCAL LO signals from any two octants that were 
separated by at least 90° in azimuth, where both octants also satisfied the lower threshold 
PRETRIGGER LO requirement. 

In addition to the high p? triggers, a prescaled sample of low bias triggers were recorded 
concurrently. These low bias triggers included beam, interaction and pretrigger events, and 
constituted ~10% of the recorded events. 



III. CHARGED TRACK RECONSTRUCTION 

The charged tracks in the selected events were reconstructed first in the PWC system. 
Since the PWC system consisted of four views, three-dimensional (3D) information was 
extracted. Space tracks were formed by combining all possible XY (UV) candidate track 
segments and searching for hits along the projections in the remaining two views. Space 
tracks were required to have a minimum of 11 hits if they involved all four PWC stations (16 
planes), 10 hits if they involved three PWC stations (12 planes), and 6 hits if they involved 
only the two upstream PWC stations (8 planes). 

The STDC pattern recognition was initially seeded by the 3D space tracks reconstructed 
in the PWC system (in order to correlate the hits in the X and Y views of the straw tubes). 
Once the correlation was performed, an iterative procedure was used to form straw tracks 
using only the hits detected in the STDC's. Straw track segments required a minimum of 4 
hits in either the XorY view. The angular resolution of straw tracks was £«0.06 mrad. After 
reconstructing the straw track segments, each space track was refit using the hit information 
from both the PWC and STDC systems. 

After identifying the downstream space tracks, track segments were reconstructed in 
the X and Y views of the SSD system. Four and five hit tracks were reconstructed, and 
then three hit tracks were formed from the previously unused hits. The SSD tracks had 
an average angular resolution that was similar to the STDC's, and an impact parameter 
resolution at the primary vertex of ~15 /zm (for p > 15 GeV/c). Track segments were 
also reconstructed from the hits in the beam SSD modules to measure the direction of the 
incident beam particle. 

Three dimensional tracks were formed in the SSD system by linking the projected 3D 
downstream (PWC and STDC) tracks to corresponding projected SSD X and Y track seg- 
ments at the center of the analysis magnet. In the bend plane of the analysis magnet, a 
link between an SSD X track segment and a downstream space track was established if the 
corrected difference between the X positions at the magnet center of the projected SSD X 
track segment and the projected downstream track was within 3.3<tax- A link between an 
SSD Y track segment and a downstream space track occurred when the corrected difference 
between the Y positions of the projected tracks at the magnet center was within 3.3<7ay and 
the corrected slope difference between the tracks was within 3.3cta6v [ 3~8"| . The matching 



resolutions for the projections (a ax and a ay) and the YZ slopes (o"A6v) were momentum 
dependent functions which were extracted from the data. SSD X and Y view track seg- 
ments were correlated with each other by virtue of being linked to the same downstream 
space track. Once the linking was complete, the primary vertex for the event was recon- 



structed [35] . Figure ^ shows a distribution of reconstructed primary vertex locations along 
the nominal beam direction [Z axis) for events accumulated during the 1990 fixed target 
run. Based upon the primary vertex location, one can identify whether the interaction oc- 
curred in the beryllium, copper, or silicon (SSD) targets. The average resolution for the Z 
location of the primary vertex was ~300 fim. After the primary vertex was located, the di- 
rection cosines, charge, and momentum of each reconstructed charged track were evaluated. 
(The momentum scale was calibrated using J/ip and Kg signals.) The average momentum 
resolution for charged tracks produced in the target region was o v jp ~ 0.0076 + 0.00026p, 
where p is the momentum measured in GeV/c. 



After the full sample of events was reconstructed and the results written in the form of 
data summary tapes (DSTs), a secondary vertex analysis was performed to search for evi- 
dence of heavy quark production. Only events that had a reconstructed primary vertex in 
the copper or beryllium targets were used in this analysis (—17.1 cm < Z < —8.5 cm). For 
each event, tracks having large transverse significance to the primary vertex were identified 
as secondary tracks. Transverse significance is defined as the measured impact parameter 
between a vertex and a given track candidate divided by the corresponding expected uncer- 
tainty. The algorithm evaluated all pairs of secondary tracks and selected only those pairs 
that were consistent with emanating from the same space point. For all such combinations, 
the algorithm determined whether any other secondary tracks had transverse significance of 
less than 3 units relative to the space point in question. All such tracks were added to the 
track list associated with that secondary vertex, and the vertex location and its associated 
uncertainties were reevaluated (via a x 2 minimization technique). A list of vertices, each 
determined by two or more space tracks was thus generated. Two track vertices were re- 
ferred to as vees and all other vertices were referred to as secondary vertices. To minimize 
the losses introduced by only utilizing secondary tracks, all space tracks in the event were 
examined to determine whether they might belong to any given secondary vertex or vee. 
This phase of the program generated a list of additional tracks that might possibly belong 
to each of the secondary vertices or vees. Neither the secondary vertices nor the vees were 
refit with any of these additional tracks. Only those events with at least one reconstructed 
secondary vertex or vee contributed to this analysis. 

IV. CHARM SIGNALS IN THE E706 DATA 

From the data sample acquired during the 1990 fixed target run, we have identified 
D°, D*^, and D^ 1 signals in fully charged modes. The D° and D*^ signals (see Fig. f|) 
do not directly contribute to the measurements presented in this report, which are based 
upon the sample of D ± mesons that were observed via their decay to the fully charged final 
state K^ir^n^. The D sample was extracted from the subset of events that contained 
at least one secondary vertex with 3 tracks or a vee with additional tracks attached to 
it. For all such vertices (vees), the three-body invariant mass was evaluated by assigning 
the charged pion mass to each of the two like-charge tracks while the oppositely charged 
track was assigned the kaon mass. To reduce the large combinatorial background, only 
secondary vertices that satisfied additional requirements contributed to the final analysis. 
The significant requirements were: (1) the impact parameter to the primary vertex of the 
parent momentum vector formed from the candidate decay products must be less than 50 /xm; 
(2) the longitudinal separation between the primary and secondary vertex normalized by the 
corresponding expected uncertainty in that separation must be at least six; (3) the impact 
parameter of each candidate decay track to the secondary vertex (vee) must be less than 
0.4 of the corresponding impact parameter to the primary vertex for every candidate decay 
track that contributed to the determination of the secondary vertex (vee) location; (4) for 
those tracks that contributed to the secondary vertex finding, the product of the secondary 
to primary vertex impact parameter ratios must be less than 0.005 for secondary vertices or 
less than 0.002 for vees. 



Figure [| shows the Kitit invariant mass spectrum from the events satisfying the above 
analysis requirements. The mass distribution contains 110 ± 17 combinations above back- 
ground in the D ± mass region (1.8 to 1.94 GeV/c 2 ) with xp = 1p z j\fs > —0.2 and 
Pt > 1 GeV/c. Figure |], which shows mass spectra in several pt intervals, illustrates 
the broad pt range populated by this D sample. This large range is a consequence of the 
E706 trigger, which preferentially selected high p? interactions. The high trigger thresholds 
coupled with the large combinatorial background compromise our ability to observe a signif- 
icant signal for D mesons with px < 1 GeV/c [as indicated in Fig. |6|(a)]. (The main source of 
background is secondary interactions in the target.) The D signal region was defined to be 
between 1.80 and 1.94 GeV/c 2 , except for the interval 1 < pt < 2 GeV/c, where a narrower 
mass range of 1.82 to 1.92 GeV/c 2 was employed. This more restrictive mass requirement 
reduces the statistical uncertainty with only a minimal loss of efficiency Since the resolution 
of the reconstructed D signal (~19 MeV) observed in the data and our MC simulation were 
consistent, the small loss of events resulting from the narrower mass range was absorbed 
into the reconstruction efficiency. The background in the signal region was estimated via a 
linear interpolation between the lower and upper sideband regions. The uncertainty in the 
number of background combinations was estimated by fitting the background to first and 
second order polynomials over several mass regions (which included the signal region). 

Figure [7] shows the tracks from an event containing a reconstructed charm particle can- 
didate with a pt of 4.1 GeV/c. The figure shows the various target elements, the first vertex 
SSD plane, and the projected charged tracks from the event which reveal a primary vertex 
and a displaced three-track secondary vertex. Due to the high transverse momentum of the 
charm candidate, the secondary vertex is well isolated from the other charged tracks in the 
event. Note that the vertical scale is magnified with respect to the horizontal scale; the 
polar angle of the widest angle track is no more than about 6° . 

V. EFFICIENCIES 

An event simulation was used to estimate the efficiency for selecting events containing 
charged D mesons and for detecting those D^ decays. The simulation includes an event 
generator and a detector simulation. The event generator simulates particle production in 
high energy collisions, and the detector simulation models the response of the detectors to 
the generated particles. The details of this simulation and the evaluation of the efficiencies 
are discussed below. 



A. Event generation 

The event generator chosen to produce full events was the pythia 5.6/jetset 7.3 pack- 
age |40| . The physics processes employed in the event generation are specified by the user. 



The physics processes investigated in this analysis included minimum bias events and a pure 
charm event sample. The former was used to tune the MC parameters to match the global 
characteristics of events observed in the data. Once the MC simulation was tuned, the 
efficiencies for triggering on and selecting events containing charm particles were evaluated 
using the pure charm event sample, pythia describes the hard scattering between hadrons 
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via leading order perturbative QCD matrix elements, and simulates the NLO contributions 
through effective K factors ||0],[|l|. Parton showers are produced via perturbative branch- 
ings of one parton into two or more partons. The pythia simulation also includes initial 
state radiation of the incoming partons, which by default, is activated, jetset handles the 
nonperturbative fragmentation of the final state colored partons into colorless hadrons using 



the LUND string model |20 ]. In addition, jetset handles the decays of unstable particles 
via a list of decay modes and branching ratios that are extracted from Particle Data Group 
tables. For each event, the event generator provides a list of the resulting stable particles 
(and their associated kinematic variables) that can be used as input to the E706 detector 



simulation. The detector simulation is based upon the geant software package [f42 |. 

The event generator was first tuned to match various distributions observed in our data 
which were relevant to this analysis. Since the trigger discriminated using electromagnetic 
depositions, it is important to reproduce the pt spectrum of particles that produce elec- 
tromagnetic showers. Figure ^] shows the spectrum of 7r°'s measured in the data and the 
corresponding spectra generated by pythia for several choices of the kj_ parameter. The 
data distribution was measured using the low bias triggers, and the MC spectrum was gen- 
erated using minimum bias events. In each case, the MC spectrum is normalized to the same 
integral as the data over the kinematic range shown in the figure. The data are reasonably 
described by the pythia result generated using a kj_ value of ~1 GeV/c, which is larger 
than the pythia default kj_ value of 0.44 GeV/c. The px spectra of charged tracks were 
also compared and exhibited a similar level of agreement. We chose a k± value of 1.0 GeV/c 
for further study. To improve the match between the charged particle multiplicity observed 
in our data and pythia, we adjusted the pythia parameter designated as the effective 



minimum transverse momentum for multiple interactions |43| ) which increased the mean 
track multiplicity by ~30%. Figure ^| shows the multiplicity distributions of charged tracks 
as reconstructed in the PWC system and in the X and Y views of the SSD system for 
both the data sample and the Monte Carlo sample (after adjustment of the aforementioned 
pythia parameter) for events that satisfied at least one of the triggers used in this analysis. 
The corresponding data and MC distributions are similar, indicating that the adjusted MC 
simulates the particle multiplicity of high energy collisions reasonably well. Rapidity distri- 
butions of charged tracks in the MC and data were also found to be in reasonable agreement. 
With these modifications, the pythia simulation adequately describes the distributions of 
final state particles observed in our data. 

B. Trigger simulation 

As previously described, the E706 trigger utilized signals from the electromagnetic section 
of the LAC. Since the EMLAC consists of ~27 radiation lengths, but only ^1 interaction 
length, photons deposit nearly all of their energy in the EMLAC, whereas hadrons usually do 
not. Consequently, high p? photons and electrons were more likely to trigger the apparatus 
than hadrons of the same pt- To investigate the impact of the requirements imposed by 
the trigger on observed events, we developed a software simulation of the online trigger. 
Corrections for the losses resulting from the high pt thresholds of the various triggers were 
evaluated by subjecting pythia events to this software trigger simulation. The main features 
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of the trigger simulation are discussed in this section. 

The trigger simulation depends upon the modeling of the energy deposited in the active 
volume of the EMLAC by incident particles. In order to generate a substantial number of 
events in a timely manner, the showers generated in the EMLAC by incident particles were 
parametrized. After accounting for the energy loss in the inactive material in front of the 
calorimeter (such as the cryostat wall), the (calibrated) energy response of the EMLAC to 
incident photons and electrons was parametrized by an energy resolution function, oe = 
0.22 2 + 0.1G 2 E + (0.01E) 2 , where E is the energy in GeV. The response of the EMLAC 



to incident hadrons was investigated using the geant software package ||42| . Individual 
hadrons at a fixed energy of 20 GeV were generated in the target region and propagated 
through the full detector simulation (resulting in full showers in the EMLAC). The simulated 
calibrated energy response relative to the incident hadron energy was evaluated for various 
incident hadron types and the resulting distributions are shown in Fig. 0(b) and (c). [For 
comparison, the corresponding distribution for incident 20 GeV photons is shown in Fig. 
|10"|(a).] These distributions were used to parametrize the response of the EMLAC to incident 
hadrons. Since these distributions were generated based upon the full shower response of the 
EMLAC, they already include the effects of the intrinsic energy resolution of the EMLAC. 
The shapes of these distributions were found to be relatively insensitive to variations in 
the incident energy for E > ^6 GeV, and hadrons below this energy were generally not 
detected. The transverse and longitudinal development of showers were also parametrized 
based upon full shower geant studies. In the transverse direction, a radial shower profile 
was used to generate the simulated energy deposition on the R strips near the centroid of 



the shower (see Fig. [11]). The longitudinal shower development for incident photons and 



hadrons was parametrized as the ratio of the energy in the front section of the EMLAC with 



respect to the total energy deposited in the EMLAC ||39| . These parameterizations of the 
response of the EMLAC to incident photons and hadrons were used to provide an efficient 
simulation of the distribution of energy deposited in the EMLAC for each of the incident 
particles generated by pythia which impinged on the active region of the EMLAC. In this 
manner, the total energy deposited on each R strip was estimated as the scalar sum of the 
energy depositions of the individual particles in that given event. 

The trigger pt detected by each R strip was evaluated by weighting the energy detected 
on that strip by the appropriate measured strip trigger weight, which increased as ~ sin^j. 
(The strip trigger weights were measured and evaluated independently for each trigger type 
used in this analysis.) From the px in the strips, LOCAL and GLOBAL px sums were 
calculated analogously to the online trigger method. These pt sums and the measured 
trigger efficiency curves for each trigger were used to evaluate the event trigger probabilities 
for each trigger type for each event. The simulated event was either accepted or rejected 
based upon these trigger probabilities. 

This trigger simulation was tested by comparing the fraction of interactions in which the 
LOCAL LO and LOCAL HI requirements were satisfied in the Monte Carlo simulation and 
the low bias data. The interactions in the low bias data were required to have a reconstructed 
vertex in the target region and events containing beam halo muons were excluded. Several 
MC event samples were generated using the minimum bias event generator and varying the 
kj_ parameter from 0.7 GeV/c to 1.3 GeV/c. For each of these samples, we measured the 
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rate at which the LOCAL LO and LOCAL HI signals were generated. Figure 0(a) shows 
the rate of LOCAL LO signal generation in the data sample compared to the corresponding 
rate determined from the MC samples. Figure |12](b) shows a similar comparison for the 
LOCAL HI signal. The open circles represent the MC results for the various kj_ choices, and 
the shaded band represents the data and its associated uncertainty. Both the LOCAL LO 
and LOCAL HI signal rates are reproduced by the pythia MC simulation using a k^ value 
of wl GeV/c. The relative rates between the SINGLE LOCAL HI, LOCAL®GLOBAL HI, 
and TWO-GAMMA triggers were also consistent between the MC and data samples. These 
observations indicate that the tuned MC simulation provides reasonable estimates of the 
rates at which high energy interactions satisfy the various triggers that contributed to this 
analysis. 

C. Trigger efficiency for charm events 

Since the tuned pythia MC simulation reproduced the kinematic distributions of final 
state hadrons, as well as the observed trigger rates in the data, we used this simulation 
to evaluate the trigger efficiency for charm events. The trigger efficiency is determined by 
calculating the probability that an event containing a D ± meson (which decays to K7nr) will 
satisfy the SINGLE LOCAL HI, LOCAL<g>GLOBAL HI, or TWO-GAMMA trigger. (We 
compared the integrated D^ cross sections determined using each trigger type individually, 
and found them to be consistent within uncertainties.) The average trigger efficiency is 
evaluated as a function of the p? of the charged D meson that decayed to the Kirir final 
state. Figure |13|(a) shows the resulting trigger efficiency for the central kj_ value of 1 GeV/c, 
as well as for other reasonable choices of kj_ values based upon Fig. [12|. In Fig. |l3|(b), 



the ratios of the trigger efficiencies for the larger (and smaller) kj_ values with respect to 
the central value are presented. The ±15% systematic uncertainty associated with this 
choice was determined based on the fractional difference between the mean value of the 
efficiencies determined via the larger and smaller kj_ choices. It is also plausible that the 
trigger efficiency will be sensitive to the value of the charm quark mass that is used in the 
MC simulation. Figure |14|(a) shows the trigger efficiency versus the transverse momentum 
of the D ± for three choices of m c (1.2, 1.35, and 1.5 GeV/c 2 ). Figure |H](b) shows the ratio 
of trigger efficiencies for the larger and smaller charm quark mass values with respect to 
the central value. The ±10% uncertainty associated with the choice of m c =1.35 GeV/c 2 
was determined based upon the fractional change in the mean value of the trigger efficiency 
when m c was varied from 1.2 GeV/c 2 to 1.5 GeV/c 2 . 

D. Tracking simulation 

In addition to correcting for losses due to the trigger, we must also evaluate the efficiency 
of reconstructing the decay vertex once the event has triggered the apparatus. All of the 
tracking detectors were modeled within the framework of geant to estimate the losses due 
to the geometrical acceptance of the tracking system and detector performance. Various 
detector effects, such as efficiency, resolution and noise, were evaluated in the data and 
then incorporated into the MC simulation. Hit efficiencies are modeled as a function of 
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transverse position in the tracking detectors. Figure [15| shows the number of hits on the 
reconstructed tracks in the PWC system, and the X and Y views of the SSD system for 
both the MC and data samples. The consistency between these distributions indicates that 
the tracking detector efficiencies are well modelled. We also tuned the MC simulation to 



reproduce the hit multiplicities in the charged tracking detectors. Figure |16| shows the MC 
and data distributions for the average hit multiplicity in the PWC planes, and the 25 /zm 
and 50 fim pitch SSD planes, respectively. Both the effects of delta rays and noise hits 



have been included in this simulation [[44]. The mean hit multiplicity from the data and 



the MC simulation are similar; however, the data distributions are slightly broader. The 
impact of increasing the hit multiplicity assumed in the MC simulation on the charged D 
meson reconstruction efficiency is discussed in the next section. Finally, we investigate how 
well the resolution of the detector is simulated. The impact parameter distribution provides 
a measure of the angular precision of the SSD tracks, and this is pertinent to separating 
heavy quark decay vertices from the primary vertex. Figure IT7] shows the MC and data 



impact parameter distributions of all charged tracks relative to the primary vertex for the 
X and Y track segments of the SSD system. The MC distribution is ~5% narrower than 
the corresponding data distribution. The effect of this difference on the estimated D^ 
reconstruction efficiency will be addressed in the next section. 

The decay of K s into 7r + 7r~ provides an opportunity to compare high statistics MC and 
data distributions relevant to secondary vertex finding. A large sample of Kg mesons were 
reconstructed in the data using displaced two track secondary vertices. A sample of full 
events enriched with K s mesons was generated using the minimum bias event generator 
from pythia, and reconstructed using the secondary vertex algorithm described in this 
paper. The background subtracted Kg signals are shown in Fig. [H|. The MC and data 
K s mass resolutions are in good agreement. Figure |T^(a) shows the 3D impact parameter 
distributions of the tc + and tc~ tracks (from the K s signal region) to the vee location. Figure 
|T5|(b) shows the difference in the reconstructed Z coordinate of the vee as determined in the 
X and Y views independently. In both cases, the MC simulation reproduces the distributions 
observed in the data. We have compared the distributions of other analysis variables as well 
and found similar agreement |3£ 



E. Charged D meson reconstruction efficiency 

The D^ reconstruction efficiency, which includes acceptance losses, detector effects, and 
analysis requirements, is evaluated using the sample of pythia charm events generated with 
k_i_=1.0 GeV/c that satisfy at least one of the high p? triggers. Figure |2UJ shows the D^ 1 



reconstruction efficiency as a function of the reconstructed D transverse momentum. The 
average reconstruction efficiency for each pt bin is defined as the value of the parametrization 



shown in Fig. |20J evaluated at the center of the bin. The reconstruction efficiency increases 
from f»9% in the lowest px bin (1 to 2 GeV/c) to ~17% in our highest pt bin (6 to 8 GeV/c). 
The inset in Fig. ^ illustrates the D reconstruction efficiency versus xf for D mesons 
with p T > 1 GeV/c. 

To estimate the systematic uncertainty in the reconstruction efficiency, two additional 
versions of the MC simulation were evaluated. One version included more noise hits than 
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the default version, and the other had reduced hit multiplicity. The extreme values were 
determined based upon variations observed in the data at different times during the data 
taking period. The variations in hit multiplicity were typically less than 10% of the central 
value. The primary causes of these variations in the average hit multiplicity are variations 
in the beam intensity and detector performance. Increasing the number of noise hits results 
in a degradation of the track and vertex resolution which reduces the D ± reconstruction 
efficiency Reducing the hit multiplicity has the opposite effect. We compared the impact 
parameter distributions of charged tracks to the primary vertex for these two MC versions 
and found that the higher hit multiplicity version more accurately reproduced the data 
result. However, the hit and track multiplicity in this version of the MC was larger than 
that observed in the data. The sensitivity of the D^ reconstruction efficiency to the choice 
of either matching the track multiplicity distribution or the impact parameter distribution 
from the MC simulation and the data reflects an uncertainty in the reported reconstruction 
efficiency. The 10% systematic uncertainty in the D ± reconstruction efficiency is based on 
the relative difference between the results from these two versions of the MC simulation. 



VI. RESULTS 

The observed D — > Ktttt signals and the efficiencies for triggering on and reconstructing 
these events are tabulated in Tables | and || for charged D mesons in the kinematic ranges 
xf > —0.2 and xp > 0.0, respectively. Using these numbers, the cross section in each, pt bin 
is given by , ) LB i where N(px) is the number of mass combinations above background in 
the given px bin, e(pr) is the efficiency for reconstructing those events, L is the integrated 
luminosity, and B is the branching ratio for D — > Kttti, which is 9.1±0.6% |45|j . The 
integrated luminosity per nucleon is 7.8±0.8 events/pb for this data sample, including target 
transverse fiducial cuts and corrections for beam absorption. In combining the signals from 
the beryllium and copper targets, nuclear effects were assumed to be negligible for the 
hadroproduction of open charm. The dependences of cross sections per nucleus on atomic 
mass are often parametrized as <JoA a , where A is the atomic mass of the target nucleus. 
Recent measurements of D meson production yield values of a consistent with one [^,46,47 



No significant pr (or xf) dependence of a is observed [0. The differential and integrated 
D ± cross sections are discussed in the following subsections. 



A. Differential cross sections 

The D ± differential cross sections per nucleon integrated over the ranges xf > —0.2 and 
xf > are presented in Table fTT| j|5| and displayed in Fig. |2T|. Due to the steeply falling 
spectra and the large widths of the pt bins, the data points are plotted at pr values, p 1 ^, that 
correspond to the average values of the cross section in the appropriate bins as determined 



from the P ythia MC pr spectra |^9[ . The uncertainty in the cross section is obtained from 
the quadrature sum of the statistical and systematic uncertainties. 

Figure |22] shows our measured D ± cross sections and the results of the pythia MC 
simulation. Common systematic uncertainties in the luminosity and branching ratio have 
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been excluded from the uncertainties shown in this figure. The shapes of the D spectra 
generated by the pythia MC simulation are consistent with our measured spectra. These 
pythia results were normalized to our measured cross section integrated over xp > —0.2 and 
1 < Pt < 8 GeV/c. The consistency between the data and the corresponding pythia result 



integrated over the range xp > [shown in Fig. 22|(b)] indicates that pythia adequately 
models the fraction of the cross section in the range xf > relative to the fraction in the 
range xf > —0.2. 



Figure ^ shows the differential charged D cross section compared to results of NLO 
pQCD calculations [pUj] . The renormalization and factorization scales employed in the cal- 
culation of these differential distributions are as follows: jjlr = /i and ftp = 2/i , where 



/io = Jrn? c + \(Pt + Pt) an d Pt and p T are the transverse momenta of the charm and anti- 
charm quarks. The solid curve shows the NLO pQCD result for charm quark production. 
The dotted curve illustrates the NLO pQCD result including the effects of Peterson et al. 
fragmentation (with e c =0.06). The other broken curves show the Peterson et al. fragmented 
NLO pQCD results for charm production supplemented with various values of (k^). The 
pQCD cross sections shown in Fig. ^3| have been normalized to our extrapolated integrated 
D ± cross section (see next section). The results are integrated over the range xp > —0.2. 
(The corresponding comparisons for results integrated over the range xp > do not ap- 
pear substantially different.) The common systematic uncertainties in the luminosity and 
branching ratio have been excluded from the uncertainties shown in this figure. The shapes 
of the theoretical pp spectra for unfragmented charm quark production are consistent with 
data on the hadroproduction of D mesons from experiments E769 [[51]] and WA82 |J . E653 
reported that the shape of the unfragmented NLO pQCD transverse momentum distribu- 
tion was "somewhat harder" than the observed distribution [|Tj|. Our data probe a larger 
Pt range, and the unfragmented theoretical pp spectrum is clearly harder than the spec- 
trum observed in our data. Introducing Peterson et al. fragmentation into the calculation 
results in a softer p T spectrum whose shape better matches the shape observed. It is clear 
that the NLO pQCD spectrum generated using m c —1.5 GeV/c 2 and supplemented with 
(k^)=3 (GeV/c) 2 is harder than our data, while the corresponding pp spectra generated 
with (k^) values between 1 and 2 (GeV/c) 2 are similar to what we observe. As noted in the 
Introduction, the D meson pp spectra reported by experiments E769 and WA82 can also be 
described by the fragmented NLO pQCD results supplemented with a (k 2 : )=2 (GeV/c) 2 [Q. 

B. Integrated cross section 

Since the shapes of the E769 and WA82 D meson pp spectra are consistent with 
the NLO pQCD calculations using Peterson et al. fragmentation supplemented with 
(k 2 : )=2 (GeV/c) 2 , we used that pQCD calculation to estimate the extrapolation factor 
necessary to account for the low pp portion (pp < 1 GeV/c) of the D^ cross section. We 
find an extrapolation factor of 2.7. The uncertainty in that factor is estimated to be 15% 
(which corresponds to the range of values obtained from the fragmented NLO calculation 
supplemented with (k 2 : ) = 1.5 (GeV/c) 2 and the unfragmented NLO pQCD calculation). An 
additional factor of 1.07±0.03 accounts for the D ± mesons produced with xp < —0.2. This 
factor was evaluated using the xp spectrum of D^ mesons from the pythia MC simulation. 
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The uncertainty was based on half the difference between the factor predicted for D + and 
D~ mesons separately. 

Using these factors, the total D ± cross section is 16.6 ± 4.5(stat) ± 4.8(syst) fib per nu- 
cleon. The systematic error estimate includes uncertainties in the trigger and reconstruction 
efficiencies, normalization, branching ratio, and the extrapolation. The D cross section for 
xp > is 11.4 ± 2.7(stat) ± 3.3(syst) fib per nucleon. This result is compared in Fig. 
24] with previous measurements P|-4T^|] of the inclusive D^ cross section in ii~ interactions 



at other beam energies. Where appropriate, the results from other experiments have been 
adjusted to reflect current branching ratio values. The error bars represent the statistical 
and systematic uncertainties added in quadrature. The figure also shows the results of NLO 
pQCD calculations of the charm cross section ||5|,|6|,|50|| . The results have been adjusted by a 
factor of 2 to account for the associated production of charm and are reduced by a factor 
of 1.6 to reflect the partial xp coverage [xp > 0) [jJ!|[52|1 . Since it is expected that the D ± 
to charm fraction is nearly constant over this energy range JI], the results are also multi- 
plied by a factor of 0.22 to account for the fragmentation process c — > D . This factor 
was extracted from published measurements of the forward production charm hadron cross 



sections ||11|| . The pQCD calculations use the SMRS2 |53| and HMRSB |54j parton dis- 
tribution functions for the pion and nucleons, respectively. The solid curves in the figure 
are the values of the calculated charm cross section generated using a charm quark mass of 
m c = 1.5 GeV/c 2 and renormalization scales of 2m c and \m c . To illustrate the sensitivity 
of the calculation to variations in the input parameters, the corresponding results are also 
displayed for a charm quark mass of 1.35 GeV/c 2 (dashed curves). The calculations are 
expected to exhibit similar sensitivities to other inputs, such as the factorization scale and 
parton distribution functions. While these theoretical uncertainties in the normalization of 
the charm cross section are large at NLO, the energy dependence of the calculated charged 
D meson cross section is less sensitive to these uncertainties and adequately describes the 
energy dependence suggested by the data points from the various experiments. 



VII. CONCLUSIONS 

We have analyzed a large sample of 515 GeV/c tt~ interactions in copper and beryllium 
targets selected via a high transverse momentum electromagnetic shower trigger to study the 
hadronic production of high transverse momentum charged D mesons. Secondary vertices 
from D^ — > ir =F 7r ± 7r ± decays were reconstructed in a charged particle tracking system which 
included silicon strip detectors, a dipole magnet, proportional wire chambers, and straw 
tube drift chambers. The data span the kinematic range xf > —0.2 and 1 < pp < 8 GeV/c, 
a range which exceeds previously reported measurements. 

The measured differential cross section is consistent with results from pythia, provided 
pythia is supplemented with a kj_^l GeV/c. Our measured charged D meson pp spectrum 
is not well described by the NLO pQCD calculations for charm quark production that 



describe the lower pp results reported by E769 |5"I |. Our data are consistent with the charm 
Pt spectra resulting from Peterson et al. fragmented NLO pQCD calculations supplemented 
with a (k^) of 1 to 2 (GeV/c) 2 . This observation is consistent with previous reports on 
the comparison of the measured inclusive pp spectra of hadroproduced charm particles to 
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fragmented NLO pQCD calculations supplemented with additional (k^) |L|. 

An extrapolation based upon our measured spectrum yields an integrated D ± cross 
section of 11.4 ± 2.7(stat) ± 3.3(syst) fib per nucleon for xp > 0.0. This value is consistent 
with the trend observed in measurements at other incident beam energies. The total D^ 1 
cross section for 515 GeV/c 7r~-nucleon interactions is 16.6 ± 4.5(sta£) ± A.8(syst) fib per 
nucleon. 
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TABLES 

TABLE I. Numbers of D^ candidates and associated efficiencies per pt bin for D^ mesons 
with xf > —0.2. For the number of D^ candidates, the first uncertainty represents the statistical 
error in the total number of combinations in the D signal region, while the second uncertainty 
was determined by varying the fitted background shapes and the range over which the background 
fits were performed. The quoted uncertainties associated with the efficiencies reflect the statistical 
and systematic errors, respectively. 



Pt bin (GeV/c) 



Number of D 



± 



Trigger Eff. (%) 



Fitted Recon. Eff. (%) 



1-2 


43 ± 12 ± 6 


0.0138 ± 0.0005 ± 0.0023 


8.9 ±0.6 ±0.9 


2-3 


43 ± 9 ± 4 


0.0559 ± 0.0028 ± 0.0093 


14.4 ±0.9 ±1.5 


3-4 


16 ±5 ±2 


0.229 ±0.018 ±0.038 


16.4 ±1.2 ±1.6 


4-6 


6±3±2 


1.49 ±0.13 ±0.25 


17.4 ±1.5 ±1.7 


6-8 


2 ±1.4 


8.81 ±2.2 ±1.5 


17.6 ±3.7 ±1.8 



TABLE II. Numbers of D^ candidates and associated efficiencies per pt bin for D^ mesons 
with xp > 0. For the number of D^ candidates, the first uncertainty represents the statistical 
error in the total number of combinations in the D^ signal region, while the second uncertainty 
was determined by varying the fitted background shapes and the range over which the background 
fits were performed. The quoted uncertainties for the efficiencies represent the statistical and 
systematic errors, respectively. 



Pt bin (GeV/c) 


Number of D^ 


Trigger Eff. (%) 


Fitted Recon. Eff. (%) 


1-2 


41 ±10 ±5 


0.0141 ± 0.0006 ± 0.0023 


11.3 ±0.8 ±1.1 


2-3 


31±7±4 


0.0581 ± 0.0033 ± 0.0096 


14.6 ±1.1 ±1.5 


3-4 


13±4±2 


0.222 ±0.020 ±0.037 


16.5 ±1.5 ±1.7 


4-6 


5±3±2 


1.36 ±0.15 ±0.23 


17.8 ±1.8 ±1.8 


6-8 





9.72 ±2.8 ±1.6 


18.5 ±4.1 ±1.9 



TABLE III. Inclusive D differential cross section per nucleon in 515 GeV/c 7r~-nucleon in- 
teractions integrated over D meson xf > —0.2 and xp > 0. The quoted uncertainties in the 
cross sections reflect the statistical and systematic errors, respectively. See the text for additional 
details. 



Pt bin 


lw 

Pt 


jg^bAGeV/c)) 


jg^b/CGeV/c)) 


(GeV/c) 


(GeV/c) 


x F > -0.2 


XF > 


1-2 


1.44±0.02 


4.9 ±1.5 ±1.2 


3.6 ±1.0 ±0.9 


2-3 


2.41±0.02 


0.75 ±0.17 ±0.19 


0.52 ±0.13 ±0.13 


3-4 


3.41±0.02 


0.060 ± 0.020 ± 0.016 


0.050 ±0.017 ±0.014 


4-6 


4.70±0.05 


0.0016 ±0.0010 ±0.0004 


0.0015 ±0.0010 ±0.0004 
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FIG. 1. Plan view of the 1990 configuration of the Meson West Spectrometer (omitting muon 
identifiers). 
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FIG. 2. Elevation view showing the configuration of the target and SSD region during the 1990 
fixed target run. Shown from left to right is the third of 3 beam SSD modules (labelled as SSD 
planes 5 and 6), the copper and beryllium targets, and the 5 vertex SSD modules (labelled as 
SSD planes 7 through 16). The instrumented regions of the SSDs are designated by the shaded 
regions. The dotted lines illustrate the size of the SSD planes. The odd numbered planes measure 
X coordinates while the even numbered planes measure Y coordinates. The strips are 50 jim wide 
on all SSD planes except for the center ±4.8 mm of SSD planes 7 and 8, where the strips are 25 ^m 
wide. Reconstructed tracks from an interaction which includes a candidate charm particle are also 
depicted in the figure. 
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FIG. 3. Longitudinal positions of the reconstructed primary vertices for a sample of events 
containing high p? showers acquired during the 1990 fixed target run. This distribution is not 
corrected for beam attenuation. 
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FIG. 4. (a) The combined K~tt + and K + ir~ invariant mass distribution for vees (two track 
secondary vertices) with pp > 1 GeV/c, and (b) the mass difference between Ktttt s and Kir 
combinations for the signal and sideband regions of the neutral D —> Kir candidates. The n s is a 
relatively low momentum (soft) pion that is attached to the primary interaction vertex. 
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FIG. 5. The K^tt tt invariant mass spectrum for those events satisfying all reconstruction 
requirements from the 1990 fixed target run data sample. All contributing Ktttt combinations have 
xf > -0.2 and pr > 1 GeV/c. 
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FIG. 7. A illustration of the reconstructed tracks in the XZ view of an event containing a 
candidate high px charged D meson decaying into the fully charged Kirir mode downstream of the 
primary vertex. 
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FIG. 8. Corrected 7r° production spectra from the E706 data and the pythia MC event gen- 
erator. The MC spectra are shown for several choices of the pythia intrinsic parton transverse 
momentum parameter kj_. 
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FIG. 10. Probability distributions of the calibrated energy response of the EMLAC for sim- 
ulated incident 20 GeV (a) photons, (b) mesons, and (c) baryons relative to the energy of the 
incident simulated particle. 
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FIG. 11. Fractional energy deposition on the R strips of the EMLAC near the shower center 
for simulated incident 20 GeV (a) photons, (b) mesons, and (c) baryons. 
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FIG. 12. The fraction of interactions in which (a) the LOCAL LO requirement was satisfied, 
and (b) the LOCAL HI requirement was satisfied, as a function of the k^ parameter of the pythia 
MC simulation. The shaded band across the plot shows the corresponding rates as measured in 
the low bias data. The width of the band is an estimate of the uncertainty in the rates measured 
in the low bias data sample. 
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FIG. 13. (a) The trigger efficiency as a function of the D pt for three values of kj_; 
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FIG. 14. (a) The trigger efficiency as a function of the D px for three values of the charm quark 
mass: m c = 1.5 GeV/c 2 , m c = 1.35 GeV/c 2 and m c = 1.2 GeV/c 2 , and (b) the ratio of trigger 
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FIG. 15. The normalized distribution of the number of hits on reconstructed tracks in (a) the 
PWC system, (b) the SSD X view, and (c) the SSD Y view. The filled circles represent the 
distribution measured in the data while the open circles are the corresponding results from the MC 
simulation. 
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(b) the vertex SSD planes that have 25 ^m regions, and (c) the other vertex SSD planes. The filled 
circles represent the distributions measured in the data and the open circles are the corresponding 
results from the MC simulation. 



37 



o 

(N 

GO 

M 
o 

S-h 

H 

O 

a 
o 

o 



0.25 
0.2 

0.15 
0.1 

0.05 




-0.03 



Data 
Pythia MC 



(a) 



_=8= 



a i i H»- 



-0.02 



-0.01 0.01 

SSD Impact Parameter - X View 



0.02 



zL 


0.3 


- 


O 




• Data 


GO 

M 

G 


0.25 


" ° Pythia MC 


rrt 






H 


0.2 


- 


<+H 








0.15 


' : 


o 












o 


0.1 


- 


en 






S-H 






Hi 








0.05 


- 



0.03 
(cm) 



(b) 



-0.03 -0.02 -0.01 0.01 0.02 

SSD Impact Parameter - Y View 



0.03 
(cm) 



FIG. 17. Impact parameter distribution of SSD tracks to the primary vertex in the (a) X and 
(b) Y views. The filled circles represent the distributions measured in the data and the open circles 
are the corresponding results from the MC simulation. 
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FIG. 18. Background subtracted ir + ir~ invariant mass distribution for vees reconstructed in the 
MC and data samples. The filled circles are the data and the open circles are the corresponding 
results from the MC simulation. The inset shows the unsubtracted 7r + 7r~ invariant mass spectrum 
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FIG. 19. (a) Spatial impact parameter distributions of the pion tracks to the Kg vee location, 
and (b) the difference in the reconstructed Z location of the vee between the X and Y views. Both 
distributions are background subtracted. The filled circles are the distributions from the data and 
the open circles are the corresponding results from the MC simulation. 
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FIG. 20. D — ► iT^vT 7T reconstruction efficiency as a function of the reconstructed D 
transverse momentum for events that satisfied at least one of the high px triggers used in this 
analysis. The curve is a parametrization of the efficiency as a function of px- The inset shows the 
D^ — ► K^tt^tt^ reconstruction efficiency as a function of xp for D^ mesons with px > 1 GeV/c 
for triggered events. Error bars reflect statistical uncertainties only. 
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FIG. 21. Differential cross section per nucleon for D production in 515 GeV/c 7r _ -nucleon 
collisions as a function of the pt of the D meson. The filled circles represent data integrated over 
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FIG. 22. Comparison of the Z? differential cross section per nucleon measured by experiment 
E706 with the corresponding result from the PYTHIA simulation MC for the kinematic ranges (a) 
xf > —0.2, and (b) xp > 0. These PYTHIA results are normalized to match our measured cross 
section integrated over xp > —0.2 and 1 < pp < 8 GeV/c. The error bars for the data represent 
statistical and systematic uncertainties added in quadrature excluding luminosity and branching 
ratio contributions. 
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FIG. 23. Comparison of the D differential cross section per nucleon measured by experiment 
E706 to the results of NLO pQCD calculations. The pure NLO pQCD unfragmented result, and the 
NLO Peterson et al. fragmented spectra supplemented with average intrinsic transverse momenta 
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FIG. 24. The integrated cross section for inclusive D production with xf > for incident 7r~ 
as a function of beam energy. Also shown are NLO pQCD calculations of the charm cross section 
in the region xp > for two choices of the charm quark mass (1.5 GeV/c 2 and 1.35 GeV/c 2 ) 
and two choices of the renormalization scale (\m c and 2m c ). In each case, the theoretical cross 
section is calculated with the factorization scale fixed at fiF = 2rn c . (Both the NA16 and NA27 
measurements were at 360 GeV, but are plotted at 356 and 364 GeV respectively, in order to 
clearly show the individual data points and their uncertainties.) 
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